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Abstract 



< 

We study bounded ancient solutions of the Navier-Stokes equa- 
tions. These are solutions with bounded velocity defined in R n x 
(—oo,0). In two space dimensions we prove that such solutions are 
either constant or of the form u(x,t) = b(t), depending on the ex- 

■ act definition of admissible solutions. The general three dimensional 
. problem seems to be out of reach of existing techniques, but partial 
| results can be obtained in the case of axi-symmetric solutions. We 

CO ' apply these results to some scenarios of potential singularity forma- 

■ tion for axi-symmetric solutions, and obtain extensions of results in a 
recent paper by Chen, Strain, Tsai and Yau [3]. 

o 



1 Introduction 

It is a well-known principle in the regularity theory of PDE that re-scaling 
procedures are very useful in studying potential singularities. For example, 
for a minimal surface £ C R n for which e E is a singular point, one should 
look at the surfaces AS in the limit A — > oo, see for example [13] . This 
"blow-up" procedure, probably first introduced by DeGiorgi in his study 
of minimal surfaces, has become indispensable in the study of singularities 
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of various geometric equations (see for example [HI [231 ES])- Analogous 
ideas were introduced in the study of many other classes of equations, such 
as semi-linear heat equations [IT] , the Navier-Stokes equations [U [7] and 
dispersive equations [HH [25], to name a few. The blow-up procedure can be 
compared to infinite magnification and therefore typically produces solutions 
of the original equation which are in some sense global. The study of such 
global solutions is often a valuable stepping stone towards understanding the 
structure of potential singularities (or the absence of singularities). In this 
paper we address some of these issues in the context of the Navier-Stokes 
equations 

u t + mVm + Vp — Au = , , 

divu = . 1 ' 

The scaling symmetry of the equations is u(x, t) — > Xu(Xx, X 2 t), p(x,t) — > 
X 2 p(Xx, \ 2 t) and can be used to "zoom in" on a solution near a potential 
singularity. There are some free parameters in this process, as we can choose 
where exactly (in space and time) we magnify (it does not have to be exactly 
at a singularity, it can for example be just before the singularity occurs), and 
which properties of the re-scaled solutions we wish to control. In this paper 
we study the situation in which we choose the L°°— norm of the re-scaled 
velocity on a certain time interval as the parameter we wish to control. The 
pressure will play no explicit role in the process. As we will see in Section El 
this leads naturally to the following global problem: 

Characterize solutions of (II. ip in R n x (— oo, 0) with (globally) bounded 
velocity u. 

Following [H], we will call solutions defined in R n x (— oo, 0) ancient 
solutions. Stated in this terminology, we are interested in ancient solutions 
of (II .ip with bounded velocity. A first guess might be that such solutions 
should be constant. To make this a plausible conjecture, one must be slightly 
more precise. Equation (II. ip has trivial non-constant solutions of the form 
u(x, t) = b(t), p(x, t) = —b'(t)x and so we need a definition of solutions which 
would eliminate these "parasitic solutions" . The right definition seems to be 
that of a mild solution (see Section [3]), which was probably introduced in 
[15j . (Implicitly it is already used in Leray's paper [21J.) Another natural 
definition often used in the literature is that of a weak solution, also essentially 
introduced in Leray's paper [21], which is defined using divergence-free test 
functions, see Section EJ This notion of solution does allow the parasitic 
solutions above. In these settings, the best possible result one can hope for 
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which is consistent with what is known about the equations would be that 
any ancient mild solution with bounded velocity is constant and any ancient 
weak solution with bounded velocity is of the form u(x, t) = b(t). We will 
prove that this is indeed the case in dimension two and also in the case of 
axi-symmetric fields in dimension three, if some additional conditions are 
satisfied (see Section E]) . The case of general three-dimensional fields is, as 
far as we know, completely open. In fact, it is open even in the steady-state 
case (u independent of t). 

The methods we use in the proofs of these results are elementary. The key 
component of the proof in dimension two is the use of the vorticity equation: 



This is a scalar equation and u satisfies the Harnack inequality (see e. g. [8]), 
which can be used to show that if u 7^ 0, then in large areas of space-time 
uo has to be almost equal to its maximum/minimum. (In fact, the strong 
maximum principle together with standard compactness results is sufficient 
to prove this.) This turns out to be incompatible with the boundedness of u. 
(One might speculate that with the condition divw = 0, a Liouville theorem 
might be true for (11.21) at a linear level, without using the relation between u 
and u. This, however, appears to be false - see [27].) The ideas behind the 
proofs of the results for axi-symmetric fields in dimension three are similar. 
In each case there is a scalar quantity satisfying a maximum principle which 
is used in a way similar to the two-dimensional case. The quantities we use 
and the corresponding maximum principles are all classical. 

There is a technical component in the proofs, since one needs to estab- 
lish that the solutions we work with have sufficient regularity. This part is 
more or less standard, and we use elementary techniques based on explicit 
representation formulae to establish the required properties. 

In the last section we use the Liouville theorems of Section [5] to obtain 
results limiting the types of singularities which may occur in axi-symmetric 
solutions of the Navier-Stokes equations. These results are inspired by the 
recent paper [1] , where significant progress in the study of the axi-symmetric 
case was made using methods quite different from the ones presented here. 
Our results on axi-symmetric singularities address some questions which were 
left open in [3] . Very recently we learned that the authors of [3] have indepen- 
dently proved results similar to those in Section using their own methods. 
Their paper (5| on the subject is expected to appear soon. 



uj t + uVu = Alv 
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It is known that axi-symmetric solutions with no swirl have to be regular, 
see [T51I29] . (We recall that the "no swirl" condition means that in cylindrical 
coordinates (r, 9, z) - see (15.51) - the ug— component of the velocity vanishes.) 
However, the case of non-zero swirl is open at the time of this writing. We 
will prove that, under natural assumptions, every potential singularity of an 
axi-symmetric solution has to be of type II, in the sense of |14| . We recall 
that a singularity of a Navier-Stokes solution u at time T is called type I if 

i / \ i 
sup \u(x,t)\ < ■ 

x ~ Vr^t 

for some C > 0. By definition, a type II singularity is any singularity which is 
not of type I. A blow up of u by a type II singularity is sometimes called slow 
blow-up, see e. g. [2]. Therefore we can rephrase our result by saying that if 
an axi-symmetric solution develops a singularity, it can only be through slow 
blow-up. We remark that Leray proved in [21] that if u develops a singularity 
at T, then 

supKx,t)| > 

for some E\ > 0. Also, the rate ^== would be the blow-up rate of a self- 
similar singularity. (It is known that these do not exist, see [22l [28] .) 

It is worth mentioning that although our results are obtained by methods 
which are more or less elementary, it seems that some of them are out of reach 
of the usual methods used in the theory of the Navier-Stokes equations, such 
as energy methods or perturbation analyses in various function spaces. This 
is because some special properties of solutions of scalar equations, although 
simple, cannot be detected at the broad level at which the usual methods 
used for Navier-Stokes are applied. A similar situation appears in the proof 
that Leray's self-similar singularities do not exist, see [221 [28], where a (non- 
classical) scalar quantity satisfying an elliptic equation is used. At the time 
of this writing, there is no known similar quantity for the general three- 
dimensional problem. 

2 Preliminaries 

Let Q C R n be a bounded domain and let T > 0. We consider the parabolic 
equation in Q x (0, T) of the form 

u t + a(x,t)Vu - Au = , (2.1) 
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with a G L£° t (fi x (0, T)). A suitable notion of a solution is for example a 
weak solution. By definition, u is a weak solution of (12.1 j) if w and V x m (the 
distributional derivative) belong to (L^ t )i oc (^ x (0, T)) and the equation is 
satisfied in distributions. It then follows from standard regularity that in 
fact ut and V 2 x u belong to {lF x t )i c(^ x (0,7")) for every p G (0, oo) and 
the equation is satisfied pointwise almost everywhere in Q x (0, T). See for 
example [20] • Therefore there is no difference between weak solutions and 
strong solutions, and we can just use the term "solution" in the context of 
(12.11) . We recall that the "parabolic boundary" of x (0, T) is <9 par (f2 x 
(0, T)) = (Clx {0}) U (dtt x [0, T}). When x G fi, the space-time points (x, T) 
belong to the "parabolic interior" of Q x (0, T) and u(x,T) is well-defined. 
We recall that the solutions of (12. ip satisfy the strong maximum principle: If 
u is a bounded solution in Q x (0, T) such that T) = sup nx ( T ) u for some 
iGll, then w is constant in Q x (0, T). In fact, a much stronger statement is 
true: non-negative solutions of (12.11) satisfy the parabolic Harnack inequality, 
see for example [8]. The Harnack inequality immediately implies the strong 
maximum principle. For our purposes in this paper the strong maximum 
principle is sufficient - we will not need the full strength of the Harnack 
inequality. Our key tool will be the following lemma which essentially says 
that the statement of the strong maximum principle is in some sense stable 
under perturbations. (This stability can be made much more precise with the 
Harnack inequality.) The lemma is certainly known in one form or another, 
but we were unable to locate in the literature the precise statement we need. 

Lemma 2.1. Let us consider equation (12. ip with bounded measurable coef- 
ficient a in Q x (0, T). Let K be a compact subset ofQ, Q' C Q' C £1 and 
t > 0. Then for each e > there exists 5 > 0, 5 = 5(Q, Q', K, T, \ \a\ t, e) 
such that if u is a bounded solution of (12. ip with sup Cx ( 0T ) \u\ = M and 
swp xeK u{x, T) > M(l - 5), then u(x, t) > M(l - e) in Q,' x (r, T) . 

Proof. We can take M = 1 without loss of generality. Assuming the state- 
ment fails for some e > 0, there must exist a sequence of coefficients a^ k \ so- 
lutions of (12.11) with a = a (fc) , and points Xk G K and (y k , t k ) G Q! x (r, T) 
such that |a (fc) | < C, < 1, u {k) (x k ,T) -> 1 and u {k) {y k ,t k ) <l-e. We 
can assume, after passing to a subsequence, that converge weakly* in 
L^° t to a, vP^ converge locally uniformly in Vl x (0, T) to G K and 

(yk,t k ) — > (y, t) G £l' x [r, T]. The regularity properties of solutions of (12.11) 
discussed above imply that u solves (12.11) with a = a, \u\ < 1 in Q x (0,T), 
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u(x,T) = 1 and u(y,t) < 1 — e. This, however, is impossible due to the 
strong maximum principle. □ 



3 Bounded solutions of the linear Stokes 
problem 

Let us first recall some basic facts about the Cauchy problem for the linear 
Stokes system, with u — (u±, . . . , u n ) : R n x (0, oo) — > R n and the right-hand 
side in divergence form: 



u t + Vp-Au = jg-fi 



dx 



Jk 



, in/Tx(0,oo) (3.1) 

divw = 

u(-,0) = M ini?" (3.2) 

Here = (fik, ■ ■ ■ , f n k) f° r k = 1, . . . ,n. Denoting by P the Helmholtz 
projection of vector fields on div-free fields and by S the solution operator 
of the heat equation, we have the well-known representation formula 



u(t) = S(t)u + S(t-s)P—f k (s)ds } (3.3) 

JO OXk 

where, as usual, u(t) denotes the function u( ■ ,t), etc. 

This can be written more concretely in terms of the kernel 

K^x,t) = (-^.A + ^-)$(x,t), 

where the "generating function" $ is defined in terms of the fundamental 
solution of the Laplace operator G and the heat kernel T by 



$(x,t)= / G(y)T(x-y,t)dy, (3.4) 

which is the same as 

See for example [T^]. Letting 

d 
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we can re-write (13.31) as 
Ui(x,t) = / T(x-y,t)u 0i (y) + / K ijk (x-y,t-s)f jk (y,s)dyds. (3.5) 

JR n JO JR n 

Note also the obvious estimates 

!^ t)|S (WW (3 ' 6) 



and 



\K ijk (x,t)\<— . (3.7) 

(l^l 2 + £) 2 



X, f 



As a consequence of (13.71) . the expression (I3.5P is well-defined for / 6 L°° 
We remark that, in contrast, solutions of 

M4 + Vp -. AM = { \ iniTx(0,oc) (3.8) 
divw = j 

u(-,0) = u iniT (3.9) 

are not well defined for / G L™ t , although the ambiguity is small. This can 
also be seen without using the explicit form of the kernel, in the following 
way: One can write, for each t, the Helmholtz decomposition of f(x,t) as 
f(x, t) = Pf(x, t) + V x 0(x, t). The projection P can be naturally defined on 
L°°(R n ) (which is mapped by P into BMO(R n ) ) only modulo constants, 
which creates an ambiguity. However, if the right-hand side is in divergence 
form, this ambiguity is cancelled by the extra derivative. 

By definition, a mild solution of the Cauchy problem (13. ip and (13.21) is a 
function u defined by the formula (13. 5p . We note that this definition does 
not involve the pressure. One can obtain (formally) an explicit formula for 
the pressure, but, unlike the formula for the velocity field u, it defines p only 
modulo a function of t (constant in x for each t) when f k is in 

The definition of mild solutions immediately implies their uniqueness. 
Also, we have standard estimates for u in terms of / = (fi, ■ ■ ■ , fn) = 
(fij)i"j=v I n particular, for u = we have the estimates 

\\u\\c^(Q(z Q ,R)) < C(a,R)\\f\\ L ~ t{R u x{0tT)) and (3.10) 
\\V x u\\ L P t{Q(z0:R)) < C(p,R)\\f\\ L ™ t{R n xi0iT)) (3.11) 
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for any a G (0,1) and p G (l,oo), where Q(z ,R) = Q((x ,t ), R) = 
B{xq,R) x (t — R 2 ,to) is any parabolic ball contained in R n x (0, T). The 
space Cp ar is defined by means of the parabolic distance y/\x — x'\ 2 + \t — t'\. 

Taking difference quotients, we see that, for u Q = 0, we have similar 
estimates for spatial derivatives: 

\\^ k M\c uiQ ( Z0> R)) < C{a,R)\\V k J\\ L ~ t{R n x{0tT)) and (3.12) 
\\^ +1 u\\ l p tiQ{zo , R)) < C(p,R)\\V k J\\ L ~ t (R n x(0,T)) (3.13) 

Moreover, a routine inspection of representation formula (13. 5 p shows that, 
when uq = 0, the time derivative satisfies, for k — 0, 1, ... , 

||V^||L^ t (i^x(0,T)) < C(T, /c)||V* +2 /||l,~ t (Jin x (o,T))- (3-14) 

We sketch the calculation leading to the last estimate in the case k = for 
the convenience of the reader: Clearly it is enough to estimate \u t (0, t)\. Let 
$ be the generating function defined in f)3.4p . which will be considered as a 
function of R n x R, with $ = for negative values of t. We can write 

Ui = (L yA $) * f jk , (3.15) 

where is a homogeneous constant coefficient operator in x of order 3 and 
* denotes space-time convolution. Applying the heat operator to (13.151) we 
can write, with a slight abuse of notation, 

(d t - A)u t = (L ijk (d t - A)$) * f jk = (L ijk G(x)5(t)) * f jk , (3.16) 

where G is the fundamental solution of the Laplacian and 5(t) is the Dirac 
distribution in t. We consider a smooth cut-off function r] = 7](x) on R n 
with j] — 1 in the unit ball B(0, 1) and r\ = outside of £>(0, 2) and set 
/' = r/f, f" = (1 — 77)/. Let us first look at u 1 ^ the contribution to U; L 
in (I3.16P coming from /'. We can move two derivatives from Lyjt to /j fe to 
obtain an estimate of (d t — A)m^(0, t) in terms of the L^ t — norm of the second 
derivatives of The estimate of (9 t — A)wf (0, t) (with the obvious meaning 
of u") is even simpler, since Lij k G is integrable in R n \ B(0, 1) and therefore 
(d t — A)u"(0, t) can be estimated in terms of the L^ t — norm of f'L. Once we 
have the estimate for (dt — A)m, the estimate for ut follows from (I3.12p . 

To define the notion of a weak solution of equation (13. ip . we follow the 
standard procedures and introduce the space Vt of smooth compactly sup- 
ported div-free vector fields (p: R n x (0,T) — > R n . We then say that a 
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bounded measurable vector field u: R n x (0, T) — > i?" is a weafc solution 
of (13 .lj) if divn = in R n x (0, T) (in the sense of distributions) and 

Jo Ir" u ^ + dx dt = Jo" Inn fk-i^V dx dt for each (p G V T . 

Lemma 3.1. For a fixed f G L™ t (R n x (0,T)) Zei u G x (0,T)) 6e 

any weak solution of ( 13.1 ft m i? n x (0,T), and denote by v the mild solution 
of the Cauchy problem ( 13. ip and ( 13.21) iwi/i no = 0. Tnen u(x,t) = v(x,t) + 
t) + where w satisfies the heat equation wt — Aw = in R n x (0, T) 
and b is a bounded measurable R n — valued function on (0, T). Moreover, we 
have the estimates 



Proof. In view of estimates ( I3.10p it is enough to consider only the case 
/ = 0. Let 0: R n x — > R be a mollifyer compactly supported in i? n x 



(-1,0), <j> e {x,t) = £- (n+1) 0(a;/£,t/£), and let u £ : R n x (0,T-e) be defined 



by n e = £ * n (space-time convolution). Let n; e be the solution of the heat 
equation in i? n x (0, T) with initial datum iy e (x, 0) = tt e (x, 0). The (smooth 
and bounded) function h e = curl(n e — w e ) satisfies the heat equation in 
R n x (0, T — e) with initial datum h e (x, 0) = and therefore it must vanish. 
Since bounded solutions of the system curlz = and divz = in R n are 
constant by Liouville's theorem, we see that u e (x,t) — w B (x,t) = b £ (t) for 
a suitable b £ : (0, T — s) — > R n . By compactness properties of families of 
bounded solutions of the heat equation we see that if e — >• along a suitable 
sequence, the functions b £ converge a. e. to an L°° function 6: (0, T) — > R n . 
The estimates follow from the constructions. □ 

Remark 3.1. In the above decomposition, the function v is of course uniquely 
determined by f , whereas the functions w and b are determined up to a con- 
stant (independent of time). In other words, the (distributional) derivative 
b'{t) is uniquely determined by u and f . 

4 Bounded solutions of Navier-Stokes 

Let us now consider the Cauchy problem for the Navier-Stokes equations: 



w\\L^ t (R«x(o,T)) < C(T)\\u\\ L ^ t ( R n x (p tT )) o,nd 

||&IU°°(0,T) < C(T)\\u\\ L ^ t ( R n x (p tT )) ■ 



(3.17) 
(3.18) 



u t + nVn + Vp — An 
divn 




(4.1) 



n(-,0) = n 



in R 



(4.2) 
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The considerations of the previous section can be repeated with f k = —u k u. 
In particular, a function u G L^ t (R n x (0, T)) is defined to be (i) a mild 
solution of the Cauchy problem (14. II) and (14. 2[) if (13.51) is valid with f k = —u k u 
and (ii) a weak solution of equation (14.11) in R n x (0,T) if divu = in 
R n x (0,T) (in the sense of distributions) and f Q f Rn u({p t + Aip)dxdt = 

Jo Jr* ~ UkU '£^^ dx dt for eacn V ^ Vr- 

Remark 4.1. obvious that the notions of weak solution and mild solution 
are also well defined under the assumption that u G L^ t (R n x (0, T") /or eac/i 
T' < T (with the possibility that ||w||x°° t (B n x(o,T')) oo as T' / Tj. 27ms zs 
a natural setting in which potential singularities of solutions of the Cauchy 
problem can be studied. Even if one considers the Cauchy problem for u 
in spaces other than L°°(R n ), such as L n {R n ) (JT^) or BMO~ l {R n ) flrfl), 
the local-in-time solution u: R n x (0,T) — > R n which is constructed for uq in 
these spaces typically belongs to L^ t (R n x (r, T — r)) for any r > 0. 

The existence and uniqueness of local-in-time mild solutions of the Cauchy 
problem (14. ip and (14.21) with u G L°° was addressed in [TU]. We briefly 
outline a slightly modified approach using standard perturbation theory. We 
define the bilinear form B : L°? t (R n x (0, T)) x L™ t (R n x (0, T)) -> L™ t (R n x 
(0,T))by 

B(u,v)i(x,t) = / -K ijk (x - y,t - s)u k (y,s)vj(y,s) dyds , (4.3) 

and we denote by U the heat extension of the initial datum it . The equation 
for w then becomes 

u = U + B(u,u) (4.4) 

and can be solved in L£° t (i2 n x (0, T)) for sufficiently small T by a fixed point 
argument, since estimate (13 ,7p easily implies 

\\B{U, v)\\ L ™ t ( R n x ( 0jT )) < CVT\\u\\Lg t (&x{Q,T))\\v\\L£ t (R"x(Q,T))- (4.5) 

We remark that (I3.10p implies that the solutions of (14.41) have enough regu- 
larity to allow us to treat (14.41) as an ODE in t, without making assumptions 
about u other than u G L™ t (R n x (0, T)). 

We recall now the regularity properties of mild solutions in L'^' t (R n x 
(0, T)). The following (optimal) result will not be needed here in its full 
generality, but we feel it is still worth mentioning: 
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Proposition 4.1. Let u E L™ t (R n x (0,T)) be a mild solution of (jPjl 
an<i (14.21) with Uq E L°° . Then for k,l = 0, 1, . . . the functions t^ +l \/^,d l t u 
are bounded and, for T' = e(k,l)\\uo\\il,r R n\ (where e(k,l) > is a small 
constant), we have 

\\$ +l V%u\\ L ~ t{R n x{or)) < C(k,l)\\u \\ L ~ iRn) . (4.6) 

Proof. This can be proved in the same way as the corresponding results 
in [12], [B] and [5J, where the authors work in function spaces other than 
L^ t . The key is an estimate of B with the same form as (14.51) but in spaces 
with norms given by the expression on the left-hand side of (14.61) . In the 
context of the L^ t (R n x (0, T))— based norms we use here, the proof is in 
fact much simpler than in that of the spaces used in the above papers, due 
to the elementary nature of estimate (14. 5p . □ 

Remark 4.2. Estimate (14. 6p says that the local-in-time smoothing properties 
of Navier-Stokes for u E L°° are the same as those of the heat equation. 
Since the solution u is constructed essentially as a power series perturbation 
around the heat extension U of u , this may not be surprising. 

Lemma 4.1. Let E L^ t (R n x (0, T)) be a sequence of mild solutions of 
(14.11) and (14. 2 j) with initial conditions wj, . Assume \ \u^\ (R n x(o,T)) < C 
with C independent of k. Then a subsequence of the sequence vS k > converges 
locally uniformly in R n x (0, T) to a mild solution u E L ( £' t (R n x (0, T)) with 

initial datum u(x,0) given by the weak* limit of a suitable subsequence of the 
(fc) 

sequence u . 

Proof. This is a routine consequence of ( 14.61) . and the decay estimate ( 13. 7[) 
for the kernel K^. □ 

We now turn to regularity properties of bounded weak solutions. Let 
u E L™ t (R n x (0,T)) be a weak solution of flHD in R n x (0,T), and let 
M = (R n x(o,T))- Let v be the mild solution of the linear Cauchy 

problem (13. ip and (I3.2p with f^ = —u^u and Uq = 0. By Lemma 13.11 we 
can write u = v + w + b with the L°°— norms of v, w and b bounded by 
N = C 1 (T)M 2 +C 2 (T)M, w t -Aw = and b is a function of t only. Hence for 
k = 0, 1, 2, . . . and 5 > the derivatives V^(u> + b) are bounded by C{k, S)N 
in R n x (5, T) by Proposition 14.11 Moreover, we have the LP— estimate (13.111) 
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for V x v. Therefore uj = curlw belongs to 1% t (Q(zo, R)) for any p E (1, oo) 
and any Q(z , R) C R n x (5,T), with 

\\ui\\ Llmz0!R)) <C(p,6,R,M). (4.7) 

Following [21], we can now use the equation for uj to obtain estimates for 
higher derivatives V x u. For n = 3 the equation for uj is 

<9 

w it - Aui = -—(ojjUi - UiUj) (4.8) 

and it is easy to check that in our situation this equation is satisfied in the 
sense of distributions. Equation (14.81) gains uj one spatial derivative in L p x t . 
The standard bootstrapping arguments and regularity estimates for harmonic 
functions now give 

\\V k x u\\L UQ ( Z0 , R) )<C(k,5,R,M) (4.9) 

for each Q(Z ,R) C R n x (8,T). Therefore, using standard imbeddings, we 
have for k = 0, 1, 2 . . . 



l|V*u||£« (Ji » x(4ir)) < C(k,S,T,M). (4.10) 

Finally, using f)3.14p we also obtain for k = 0, 1, 2 . . . 

\\V%(u-b)\\ L ^ t{R n xis>T)) <C(k,S,R,M). (4.11) 

(We adopt the usual convention that the value of C can change from line to 
line.) 



5 Liouville theorems 

Let us first consider the Navier-Stokes equations in two space dimensions. 

Theorem 5.1. Let u be a bounded weak solution of the Navier-Stokes equa- 
tions in R 2 x (— oo, 0). Then u(x, t) = b(t) for a suitable bounded measurable 
b: (-oo,0) -> R 2 . 
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Proof. In two space dimensions the vorticity is a scalar quantity defined by 



U = U 2 ,l - Ml,2 , (5.1) 

where the indices after comma mean derivatives, i. e. 1*2,1 = ~^ u 2, etc. By 
the results of Section HI the function ui is uniformly bounded together with its 
spatial derivatives. Moreover, its time derivative is also uniformly bounded. 
The vorticity equation in dimension two is 

u t + uVlu - Aw = . (5.2) 

Let Mi = sup H 2 x (_ oo ) u, M 2 = inf R 2 x (_ oo ) u and assume that Mi > 0. By 
Lemma [2J] there exist arbitrarily large balls Qr = Q((x,i), R) = B(x,R) x 
(?— R 2 ,i) such that uj > M\/2 in Q((x, t), R). For such balls we have 

/ ujdxdt>-nM x R^. (5.3) 

JQr 

On the other hand, denoting by n the normal to the boundary of B(x,R), 
we can also write 

J QR ujdxdt = f QR (u2,i-u 1>2 )dxdt = 

JdB(x,R)x(t-R2,t)( U 2 n l ~ M l™2) dsdt<CR\ 

Clearly (15.31) is not compatible with (15.41) . unless Mi < 0. In the same way we 
conclude that M2 > and therefore u must vanish identically. Hence curl u = 
in R 2 x (—00, 0) which, together with divw = and the boundedness of u, 
implies (by the classical Liouville theorem for harmonic functions) that u is 
constant in x for each t. □ 

It is not known if a result similar to Theorem 15.11 remains true in three 
spatial dimensions. In fact, the problem is open even in the steady-state 
case. However, under the additional assumption that the solutions are axi- 
symmetric, one can obtain some results which seem to be of interest. We 
recall that a vector field u in i? 3 is axi-symmetric if it is invariant under rota- 
tions about a suitable axis, which is often identified with the £3— coordinate 
axis. In other words, a field u is axi-symmetric if u(Rx) = Ru(x) for every 
rotation R of the form 



R 
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In cylindrical coordinates (r, 9, z) given by 



Xi = rcos9, x 2 = rsm9, x 3 = z, (5.5) 

the axi-symmetric fields are given by u = u r j^ + + u z j^, where the 

coordinate functions u r , ue, u z depend only on r and z. In these coordinates, 
the Navier-Stokes equations become 

U rt + U r U r r + U z U r z h p r = L\U r (5.DJ 

, UrUe A Ud fK ^ 

u et + u r u e , r + u z uq iZ H = Auo n (.5.7; 

u zt + u r u Z} r +u z u ZjZ +p jZ = Au z (5.8) 
(ru r ) r 

^^ + u z , z = 0, (5.9) 



where A is the scalar Laplacian (expressed in the coordinates (r, 9,z)), u 



i "r. z 



denotes the partial derivative §^u r , etc. The equation for ug is of special 
interest, as it is decoupled from the pressure. The role of the non-linear 
terms in this equation can be seen by considering the inviscid case (Euler's 
equations), wherein equation (15.71) is replaced by 

U r Ug 

u et + u r u d ,r + u z ue, z -\ = 0, (5.10) 

r 

which is the same as 



{ru e )t + u r (ru e ),r + u z {ru e ), z = 0. (5.11) 

Equation (15.111) says that the quantity tuq "moves with the flow". This is 
a special case of Kelvin's law that the integral of Uidxj along curves moving 
with the flow is constant. In the situation considered here, the curves are 
circles centered at the x 3 — axis and lying in planes perpendicular to it. 
In view of (15.111) . it is natural to re- write (15.71) as an equation for rug: 

2 

{ru e ),t + u r (rug), r + u z (ru e ) tZ = A(ru ) - -(rue), r (5-12) 

The infinitesimal version of Kelvin's law, which is Helmholtz's law that 
vorticity "moves with the flow" (for inviscid flows), gives in the case of axi- 
symmetric flows without swirl (uq = 0) another quantity which moves with 
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the flow, namely -^ £ . Here u = curlw, as usual, and in cylindrical coordinates 
we write u = u r j^ + ug^ + uj z j^. (For axi-symmetric flows without swirl we 
have uj r = 0, uj z = 0, and we can write uj = ^e^g- Therefore the situation is 
similar to two-dimensional flows.) 

Hence for axi-symmetric solutions of Euler's equations without swirl we 
have 

(^)* + «r(^),r + «.(^)„ = 0. (5.13) 

ry f f 

This is nothing but the 6 1 — component of the equation for u, and can be 
of course obtained by simple calculation, without any consideration of the 
Helmholtz law. For axi-symmetric solutions of Navier-Stokes without swirl 
the last equation becomes 

fVg. ,UJg UJg UJg 2 UJg . 

( — ) t + u r ( — ) r + u z ( — ) z = A( — ) + -( — ) r . (5.14) 

ty <y <y fy ry <y> 

Remark 5.1. For a smooth vector field u, the apparent singularity of 
is only an artifact of the co-ordinate choice. The quantity — is actually a 
smooth function, even across the x^ — axis, as long as u is smooth. 

The diffusion term on the right-hand side of equation ( 15. 14ft can be inter- 
preted as the 5— dimensional Laplacian acting on 5*0(4)— invariant functions 
in R 5 . We write r = \/y\ + y\ + y\ + yj, y<s = z and we note that for 
f(yi, ■ ■ ■ , Vh) = f(r, z) we have 

^f(K,...,V*) = (^ + ^+0±){r,z). (5.15) 

Therefore, with a slight abuse of notation, we can write the equation (15.141) 

as 

(^) t + u r (^), r + u z (^), z = A 5 (^). (5.16) 

fy If f f 

Theorem 5.2. Let u be a bounded weak solution of the Navier-Stokes equa- 
tions in R 3 x (— oo, 0). Assume that u is axi-symmetric with no swirl. Then 
u(x,t) = (0, 0,&3(i)) for some bounded measurable function 63: (— 00, 0) — > 
R. 

Proof. The idea of the proof is the same as in the two-dimensional case. 
By the results of Section HI we have |Vjjhz| < Ck in R 3 x (— 00, 0), and this 
implies that ^ is bounded in R 3 x (— 00, 0). Let Mi = sup J?3x („ OOi0 )(-7 a ) 
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and assume Mi > 0. Applying Lemma I2TT1 to equation (15.161) . considered 
as an equation in R 5 x (— oo,0), we see that > M\j "2 in arbitrarily large 
parabolic balls (with suitably chosen centers). However, this would mean 
that ujg is unbounded, a contradiction. Therefore M\ < 0. In the same way 
we show that M 2 = inf fi 3 X (_ oo ) — > 0, and hence ug vanishes identically. 
For axi-symmetric vector fields with no swirl this means that to = and 
the proof is finished by again applying the Liouville theorem to the system 
curlw = 0, divw = . □ 

The validity of Theorem 15.21 in the absence of the "no swirl" assumption 
is still an open problem. The following theorem, however, is a partial result 
in that direction: 

Theorem 5.3. Let u be a bounded weak solution of the Navier- Stokes equa- 
tions in R 3 x (—oo,0). Assume that u is axi-symmetric and, in addition, 
satisfies 

C 

\u(x,t)\ < in R 3 x (-oo,0). (5.17) 

v x i + x l 

Then u = in R 3 x (— oo, 0) . 

Proof. We will use the cylindrical coordinates (r, 8, z) given by (15.51) . We 
set / = rug and recall that 

ft + U r f,r + U z f, r = Af - -f jr . (5.18) 

r 

For A > we let f x {x,t) = f(Xx,X 2 t) and u x (x,t) = Xu(Xx,X 2 t). We note 
that f x again satisfies (I5.18P with u replaced by u x , a consequence of the fact 
that u x satisfies Navier-Stokes. Under our assumptions we have 

|/ A | < C in R 3 x (-oo, 0) uniformly in A > 0, and (5.19) 
C 

\u x \ < — in R 3 x (-oo, 0) uniformly in A > 0. (5.20) 

r 

Let M = sup i?3x (_ 00 o) /. We will show that M < 0. Arguing by contradic- 
tion, let us assume that M > 0. Let us fix some 5 > 0. (It is instructive to 
think of 5 as being small, although one can also take 5 — 1, for example.) 
By re-scaling / — > f x we can move points where f x is "almost equal to M" 
close to the x%— axis. Using this and Lemma 12.11 we see that for any (large) 
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T\ > 0, L > and R > and any (small) e > we can find A > such that 
f x > M — e in a space-time region of the form 



TZ± = {i£ i? 3 ,5 < \J x\ + x\< R, —L + x 3 < x 3 < L + x 3 } x (t — Ti,i) . 

(5.21) 

Consider a smooth axi-symmetric cut-off function tp(x,t) supported in 



{x G i? 3 , yxf + < R, -L + x 3 < x 3 < L + x 3 } x (t-T u t) (5.22) 
such that if = 1 in 



{x G i? 3 , + < i?-l,-L + l + x 3 <x 3 < L-l + x 3 }x(t-Ti + l,t-l) 

(5.23) 

and, moreover, \ip t \ < 1, \(p }T \ < 1 and |(/?^| < 1 everywhere. (A natural 
choice is, for example, <p(r,z,t) = £,(r)rj(z)((t) for suitable functions £, r\, ( 
of one variable.) Multiplying the equation for f x by and integrating over 
space-time, we obtain 

f / (f t X + u\f x +u x z f x -Af x )<pdxdt=[° [ --f x r! fdxdt. 

J-oo JR3 J-oo Jr? r 

(5.24) 

This equality will be shown to be impossible when M > 0. In the integral 
on the left-hand side of (15. 24ft one can change f x to f x — M and integrate 
by parts to obtain 

I [ (f x -M)(-^ t -u x V^-A v )dxdt = I + II + III. (5.25) 

J-oo JR' A 

We have f x — M = 0(e) in 7Z\ and therefore, if we allow correction terms 
of size 0(e), we can restrict the spatial integration in these integrals to the 
region {yx\ + x\ < 5}. Using (I5.19P and (15.201) . it is not hard to see that 

|/| < CL5 2 + 0(e), \II\ < C5T X + 0(e) and \III\ < C5 2 T X + 0(e). 

(5.26) 

(We remind the reader that the value of C can change from one expression 
to another.) On the other hand, the right-hand side of (I5.24p can be written 
as follows: 



/0 r 2 /*0 roo poo 

/ — f x tpdxdt = An l / / f X Lp, r drdzdt 
■oo JR 3 r ' J-oo J-oo Jo 



(5.27) 
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The key point then is that f x vanishes at the x 3 — axis and is equal to M+0(e) 
on most of the support of (p >r . It is easy to check that the last integral in 
( I5.27P is equal to 

/oo /»oo 
/ <p(0,0,x 3 ,t)dx 3 dt + O(e) < -torMLTi + C^ + CL+Oie) . 
-00 J —00 

(5.28) 

For M > 0, this leads to a contradiction to (15.261) and (I5.24p if L and T\ are 
sufficiently large and e is sufficiently small. We have proved that sup / < 0. 
It follows in a similar way that inf / > and therefore / must vanish. This 
means that the solution u is swirl-free and we can apply Theorem 15.21 to 
conclude that u = in R 3 x (—00, 0). □ 



6 Singularities and ancient solutions 

We will now consider the consequences of an assumption that a singularity 
exists in a solution of the Cauchy problem for the Navier-Stokes equations 
(14.11) and (14.21) . We aim to show that singularities generate bounded ancient 
solutions, which are solutions defined in R n x (— oo,0). More precisely, an 
ancient weak solution of the Navier-Stokes equations is a weak solution de- 
fined in R n x (—00, 0), and u is an ancient mild solution if there is a sequence 
T\ — > —00 such that u{ ■ , TJ) is well-defined and u is a mild solution of the 
Cauchy problem in R n x (TJ, 0) with initial datum u{ ■ ,Tf). (We remark that 
even if u is a bounded weak solution of Navier-Stokes in R n x (—00, 0), the 
function u(- ,t) may not be well-defined for each t, see Section HI On the other 
hand u{ ■ , t) is well defined for almost every t for any u G L^ t (R n x (—00, 0)).) 

Lemma 6.1. Assume that Ui, I — 1,2,... is a sequence of bounded mild 
solutions of Navier-Stokes defined in R n x (TJ, 0) (for some initial data) with 
a uniform bound \ui\ < C, andT\ \ —00. Then we can choose a subsequence 
such that along the subsequence the ui converge locally uniformly in R n x 
(—oo,0) to an ancient mild solution u satisfying \u\ < C in R n x (— 00, 0). 

Proof. This is an easy consequence of the results in Section HI □ 

Remark 6.1. Another easy result, which is nevertheless a useful addendum 
to the Liouville theorems of Section^ is the following: A bounded ancient mild 
solution u(x, t) of the Navier-Stokes equations which is of the form u(x, t) = 
b(t) is constant (independent oft). 



18 



We leave the proof of the last statement to the reader as a simple exercise. 

Recall from Section @] that for any u G L°°(R n ) the Cauchy problem (14. ip . 
(14.21) has a unique local-in-time mild solution u. Assume now that the mild 
solution develops a singularity in finite time, and that (0,T) is its maximal 
time interval of existence. Let h{t) = sup xeRn \u(x,t)\. By a classical result 
of Leray ([21]) we have 

m > ^= ("I 

for some E\ > 0. Let H (t) = sup 0<s<4 h(s). It is easy to see that there exists 
a sequence t k /* T such that h{t k ) = H{t^). Let us choose a sequence of 
numbers ^ k \ 1. Let N k = H(t k ) and choose x k € R n such that M k = 
\u(x k ,t k )\ > N k /j k . Let us set 



^..J-^ + i.fc + Jj). (6.2) 



The functions are defined in R n x (A k ,B k ), with A k = —M^t k and 
B k = Ml{T - t k ) > ej-?l and satisfy 

\v^\ < 7fc in R n x (A k , 0) and \v {k) (0, 0) | = 1. (6.3) 

Also, are mild solutions of the Navier-Stokes equations in R n x (A k , 0) 
with initial data v^(y) = j^-u (x k + -^-). By Lemma [6TT| there is a subse- 
quence of converging to an ancient mild solution v of the Navier-Stokes 
equations. By our construction, we have \v\ < 1 in R n x (— 00, 0) and 
l^(0,0)| = l. ' 

We have proved the following statement: 

Proposition 6.1. A finite-time singularity arising from a mild solution gen- 
erates a bounded ancient mild solution which is not identically zero. 

Without further information about the situation at hand, the proposi- 
tion may not be very useful. By itself, the existence of non-zero bounded 
ancient solutions is not surprising. (Consider constants, for example.) How- 
ever, if (non-zero) constant solutions can be excluded (for example by a 
scale-invariant estimate) and a Liouville-type theorem for ancient solutions 
is available, then finite-time singularities can be ruled out. 

A simple example of such a situation arises in the context of the Ladyzhen- 
skaya-Prodi-Serrin regularity criterion. Assume that a finite T > is the 
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maximal time of existence of a mild solution (with a suitable initial condi- 
tion). Let p, q G (l,oo) with n/p + 2/q = 1. Then | \u\ |l p '«(r™x(o,t)) = +oo- 
To see this, it is enough to note that if the L^'f— norm of u was finite, the 
function v constructed by the above procedure would have to vanish identi- 
cally a. e. due to the invariance of the L^' q t — norm under the scaling used 
in the procedure, along with the fact that the finiteness of the L^'f— norm 
implies its "local smallness". But v has to be smooth (by the results of 
Section HJ) and |u(0, 0)| = 1, a contradiction. 

A more interesting application of the procedure gives Theorems 16.11 and 
16.21 below, which can be thought of as generalizations of recent results in [4] . 

Theorem 6.1. Let u be an axi- symmetric vector field in R 3 x (0, T) which 
belongs to L^ t (R 3 x (0, T')) for each T' < T. Assume that u is a weak 
solution of the Navier-Stokes equations in R 3 x (0, T) and that 

\u(x,t)\ < , c o -- mR 3 x (0,T). (6.4) 

y/xj+x* 

Then \u\ < M = M(C) in R 3 x (0, T). Moreover, u is a mild solution of the 
Navier-Stokes equations (for a suitable initial datum). 

Remark 6.2. By the results of Section^ regarding mild solutions we see that 
u is in fact smooth in R 3 x (0, T) with pointwise bounds on all derivatives in 
R 3 x (r, T) for any fixed r > 0. 

Proof. We first prove the statement assuming that u is a mild solution (for 
a suitable initial datum). This situation is in fact the main point of the 
theorem. The fact that we can weaken the assumptions from mild solutions 
to weak solutions in the formulation of the theorem (while keeping the other 
assumptions the same) is only of marginal interest. 

Arguing by contradiction, let us assume that u is a mild solution which 
is bounded in R 3 x (0, T') for each T' < T and develops a singularity at 
time T. We now use the re-scaling procedure described in the paragraph 
preceding Proposition 16.11 to construct a bounded ancient mild solution v. 
Let Xk and M k be as in the construction. We will write Xk = (x' k ,X3k), 
with x' k = (x\k,X2k)- An obvious consequence of assumption 16.41 is that 
\%'k\ — ~w- This implies that the functions v^ k \y,s) are axi-symmetric with 
respect to an axis parallel to the y%— axis and at distance at most C from it. 
Therefore we can assume (by passing to a suitable subsequence first) that the 
limit function v is axi-symmetric with respect to a suitable axis. Moreover, 
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since assumption (16.41) is scale- invariant, it will again be satisfied (in suitable 
coordinates) by v. Applying Theorem 15.31 and using (16.41) we see that v = 0. 
On the other hand, \v(0, 0)| = 1, a contradiction. This finishes the main part 
of the proof. 

It remains to show that, under the assumptions of the theorem, u is a 
mild solution. To do this we inspect the decomposition of u constructed in 
Lemma [3. II with fa = —u^u. Using the decay of the kernel (13. 7p and of the 
heat kernel, it is easy to check that, under the assumption (16.41) . all the terms 
in the decomposition u = v + w + b will again satisfy (16. 4p . It follows easily 
that b must vanish and therefore a is a mild solution. □ 

Theorem 16. II can be used to prove the following result: 

Theorem 6.2. Let u be an axi- symmetric vector field in R 3 x (0,T) which 
belongs to L^ t (R 3 x (0,T')) for each X" < T. Assume that u is a weak 
solution of the Navier-Stokes equations in R 3 x (0, T) satisfying 

\u\ < ^£= f in R 3 x (0,T). (6.5) 

In addition, assume that there exists some Rq > such that 

for \ Aq + x\ > Ro and < t < T we have \u(x, t)\ < , = , (6.6) 

as is for example the case when u is a mild solution with initial datum uq 
decaying sufficiently fast at oo. 

Then \u\ < M — M(C) in R 3 x (0,T). Moreover, u is a mild solution of 
the Navier-Stokes equations (for a suitable initial datum). 

We remark that the statement fails, for trivial reasons, if we drop assump- 
tion (I6.6p . (Consider u(x,t) = b(t).) The fact that (16.61) is satisfied when uq 
decays sufficiently fast at oo (e. g. when it is compactly supported) follows 
for example from [Tj [2] . 

Proof. We have seen in the proof of Theorem 16.11 that (16. 6p implies that u is 
a mild solution for a suitable initial datum and is therefore smooth in open 
subsets of R 3 x (0,T). We define 

f(x,t) = \x'\ \u(x,t)\ = J~x\ + x\ \u(x, t)\, (6.7) 

where, as above, x' = (xi,^)- By Theorem 16.11 it is enough to prove that 
/ is bounded in R 3 x (0,T). Let h{t) = sup R3 f(x,t), H(t) = sup 0<T<t h{r). 
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Assume / is not bounded and choose t k /* T and G i? 3 such that = 
f{xk,t k ) = h(t k ) = H(t k ) / oo. Let \ k = \x' k \ and, for y G -R 3 , s G 
(-TA^ 2 ,0), define 

t; (fc) (y, S ) = i; (fc) (l/ , , 2 /3,s) = A feM (A fc y , ,A fc2 /3 + a;3 fc ,T + A^) . (6.8) 

We note that the sequence A& is bounded due to (16. 6p . Set = —(T — tk)\ k 2 . 

Since ( 16 ,5p is invariant under the Navier-Stokes scaling, the functions v ^ 
satisfy 

{k) \ < -S= m R 3 x (-TA fc 2 , 0) , (6.9) 



\v 

where C is the same as in (16.51) . 

Moreover, from the construction we have 

\v {k) (y,s)\<^ mR 3 x(-TX^,s k ). (6.10) 

Note also that by the elementary inequality min(l/a, 1/6) < 2/(a + 6) , 
estimates (16.91) and (I6.10p imply 

\vO\v,s)\< *L M *. inR 3 x(-T\-\s k ). (6.11) 
M fcv /-s + C\y'\ 

Let 7 C i? 3 be the unit circle {y G -R 3 , \y'\ = 1,2/3 = 0}. We have, by 
construction, \v^( • ,s k )\\ = M k which, together with ( 16.91) shows that s k > 

-C 2 M~ 2 . 

Therefore, roughly speaking, as k — > 00, the sequence v^ k > blows up along 
7. If we knew that the satisfied local energy estimates with bounds 
independent of k, the blow-up along 7 would be in contradiction with the 
partial regularity theory in [3], since the one-dimensional Hausdorff measure 
of the blow-up set must be zero. One can in fact work along these lines and 
finish the proof, but the procedure is not simple. 

One can alternatively finish the proof by another scaling argument (one 
could do both scalings in one step, but the two-step procedure seems to 
be more transparent): Denoting by e% the vector (1,0,0), for x G R 3 and 
r G (A k , 0] where A k = M|(— TA^ 2 - s k ), we define 



„w (l , T) = ^« (ei + JL 4t + _L ) . (6 . 12) 
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We will consider the cylinders 



C k = {x G R 3 , yj(x x + M k f + xl<^}. (6.13) 



It follows from our definitions that 



\w 



(k) 



(0,0)| = 1 and \w {k \x,r)\ < 2 in (R 3 \C k ) x (A k ,0) . (6.14) 



Note also that (16.111) implies 



\w ik) (x,r)\ < k = in C k x (A k ,0) (6.15) 

M - M ky /-f + Cv /(x 1 + M k y + xl 

and that (16.91) implies 

|w (fc) (a;,r)| < —= in R 3 x (A k , 0). (6.16) 

Since the functions are mild solutions of the Navier-Stokes equations in 
(A k ,0) (for suitable re-scalings of the initial datum w ), in view of bound 
(16. 161) we can choose a subsequence of the sequence w^ k \ which we again 
denote by w^ k \ such that the w^ k ' converge uniformly on compact subsets 
of R 3 x (— oo, 0) to an ancient mild solution w. In view of (I6.14p we have 
\w\ < 2 in R 3 x (— oo, 0). Moreover, since the solutions v ^ are axi-symmetric 
and M k /* oo, it is easy to see that w is independent of the x<i~ variable. 
Applying Theorem 15.11 and Remark 16.1 1 to the field (wi,w 3 ), we conclude 
that (wi,W3) must vanish identically, and this easily implies that w = in 
R 3 x (—oo,0). This would give a contradiction with \w^ k '(0, 0)| = 1 if we 
could prove that w^ k \0, 0) — > w(0, 0), which is not immediately obvious since 
our bound of sup^, \w^ k '(x, t)\ may not be uniform as r — > 0. However, by 
( I6.14p the only possible problem may occur due to the contribution from the 
cylinder C k . In the cylinder we can use the bound (I6.15P to show that the 
contribution of the dangerous part of to the representation formula (13.51) 
is negligible (in the limit k — > oo). Applying the representation formula (13.51) 
in R 3 x (—1,0) with w^ k \x, — 1) as initial datum and fji = — iUj iuj and 
using the bound (I6.15P together with the decay of the kernel (13. 7p . one sees 
that it is enough to estimate the integral 

J(M)= f I .. ; r-\ Wh2 (M1 [ ,, 2 dx'dx 3 dr. (6.17) 
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An easy calculation shows that I(M) — > as M — > oo. This shows that the 
contribution from the region where \w^\ > 2 to the representation formula 
(13.51) (with fji = —WfiUj) is negligible (in the limit k — > oo) and therefore 
(by (13.10P ) the sequence converges to w uniformly in B(0, 1) x [—1,0]. 
Therefore \w(0, 0)| = 1, which gives the sought-after contradiction. 

□ 
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